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Abstract 

All quasilinear NLS-type systems with fourth order symmetry are listed. Vector 
generalizations of some of them are constructed. Local master symmetries for several 
systems from the list are found. 

1 Classification result. 

It is well known that the following class of systems of evolution equations 

u t = u xx + F(u,v,u x ,v x ), , . 

v t = -v xx + G(u, v, u x , v x ), 

is very rich in integrable cases. In the papers [l]-[5] by Mikhailov, Shabat and Yamilov, 
all systems (1), possessing higher conservation laws, were classified. Hence, these authors 
have found all systems (1) that can be integrated by the inverse scattering method (S- 
integrable equations in the terminology by F. Calogero [6]). 

However, there are integrable cases that are not in their classification. As an example, 
consider the system 

Ut — u xx — 2uu x — 2vu x — 2uv x + 2u 2 v + 2uv 2 , ,„x 
v t = —v xx + 2vu x + 2uv x + 2vv x — 2u 2 v — 2uv 2 , 
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first discussed in [7] (see [8] for generalizations). It can be reduced to 



U t = U i 



v t = -v x 



XX 



by the following substitution of the Cole-Hopf type 



U x 



(3) 



u = 



(u + vy 



V = 



(u + vy 



The system (2) has no higher order conservation laws, but it has higher order symmetries. 
This is a typical feature of linearizable systems like the Burgers equation (C-integrable 
equations). Therefore, it would be interesting to classify the systems (1), which have 
higher order symmetries. As a result, all S-integrable and C-integrable systems would be 
found. 

The complete classification problem is very difficult. Here we consider only the most 
interesting (from our opinion) subclass of systems (1). Namely, we consider equations 
linear in all derivatives of the form 

(u t = u xx + Ai(u, v)u x + A 2 (u, v)v x + Aq(u, v) ,^ 
1 v t = -v xx + Bi(u, v)v x + B 2 (u, v)u x + B Q (u, v). 

without any restrictions on the functions Ai(u, v), Bi(u, v). We apply to such systems the 
simplest version of the symmetry test (see [9]- [12]). 

Lemma. If system (4) has a fourth order symmetry 



where A and -Bo are polynomials of at most fifth degree. 

The coefficients of the last system satisfy an overdetermined system of algebraic equa- 
tions. The most essential equations are 




(5) 



then the system is of the following form 



u t = u xx + (a 12 uv + aiu + a 2 v + a )u x + (p 2 v + p n u 2 + piu + p )v x + A (u, v), 
v t = ~v xx + (b 12 uv + biv + b 2 u + b )v x + (q 2 u + q u v 2 + q t v + q )u x + B (u, v), 



P2(b 12 - q n ) = 0, 
^2(^12 - qn) = 0, 



^2(^12 -P11) = 0, 
q 2 (a 12 -pu) = 0, 



p 2 (a 12 + 2b 12 ) = 0, 
q 2 (b 12 + 2a 12 ) = 0. 



012(012 - 012 + qn - Pu) = 0, 
(012 -Pn)(pn - 9n) = 0, 
(a 12 -p u )(a 12 - b 12 ) = 0, 



012(012 - 612 + Qn - Pu) = 0, 
(bi2 ~ qu)ipu - qn) = 0, 
(612 - 5n)(ai2 - 612) = 0. 
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As usual, such factorized equations lead to a tree of variants, which was investigated by 
the computer algebra program CRACK 1 [19]. 

Solving the overdetermined system, we don't consider so called triangular systems like 
the following 

v* = u xx + 2uv x , v t = -v xx - 2vv x . (6) 

Here the second equation is separated and the first is linear with the variable coefficients 
defined by a given solution of the second equation. 

Theorem. Any nonlinear nontriangular system (4), having a symmetry (5), up to scal- 
ings of t, x, u, v, shifts of u and v, and the involution 

u <-> v, t <-> —t (7) 

belongs to the following list: 

( u t = u xx + (u + v)u x + uv x , , g , 
\v t = -v xx + (u + v)v x + vu x , 



ut = u xx — 2{u + v)u x — 2uv x + 2u 2 v + 2uv 2 + au + f3v + 7, 
v t = —v xx + 2{u + v)v x + 2vu x — 2u 2 v — 2uv 2 — au — f3v — 7, 



(9) 



u t = u xx + vu x + uv x , 
v t = -v xx + vv x + u x , 



(10) 



Ut = u xx + 2vu x + 2uv x + 2uv 2 + u 2 + au + f3v + 7, 
. v t = -v xx - 2vv x - u x , 



u t = u xx + av x + (u + v) 2 + (3(u + v) + 7, 
v t = -v xx + au x - (u + v) 2 - (3(u + v) - 7, 



u t = u xx + (u + v)u x + Aav x + a(u + v) 2 + j3(u + v) + 7, 
ft = -fxx + (u + v)v x + 4au x - aiu + v) 2 - (3(u + v) - 7, 



M t = ^xx + 2cra 2 t> :E + 2(5uvu x + ct(/3 — 2a)-u 3 t> 2 + 7-u 2 f + 5u, 
v t = —v xx + 2av 2 u x + 2(5uvv x — a(f3 — 2a)u 2 v 3 — ^uv 2 — 5v, 



(12) 



(13) 



(14) 



1 Wolf, T. and Brand, A., The Computer Algebra Package CRACK for Investigating PDEs, manual 
+ software in the REDUCE network library (1992) 
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u t = u xx + 2uvu x + (a + u 2 )v x , , , 

v t = -Vxx + 2uvv x + ((3 + v 2 )u x , 



( u t = u xx + 2auvu x + 2au 2 v x — af3u 3 v 2 + ju, 
I v t — —v xx + 2(3v 2 u x + 2f3uvv x + af3u 2 v 3 — 71;, 

u t = u xx + 2uvu x + 2(a + u 2 )v x + u 3 v 2 + f3u 3 + auv 2 + ju, 
v t = —v xx — 2uvv x — 2{(5 + v 2 )u x — u 2 v 3 — f3u 2 v — av 3 — jv, 



u t = + 4uvu x + 4:U 2 v x + 3vv x + 2u 3 v 2 + uv 3 + au, 
v t = —v xx — 2v 2 u x — 2uvv x — 2u 2 v 3 — v 4 — av, 

f u t = + 4mm x + 2m> x , 

\ f t = — f^i: — 2t> u s — 2-uf x — 3u 2 v — v 3 + av, 

(u t = u xx + vv x , 
I Vt = —v xx + U X) 



u t = u xx + 6(u + v)v x - 6(u + v) 3 - a(u + v) 2 - (3(u + v) — 7, 
Vt = ~v xx + 6(u + v)u x + 6(u + v) 3 + a(u + v) 2 + f3(u + v) + 7, 



V* = -^xx + UU X . 



(16) 

(17) 

(18) 
(19) 
(20) 

(21) 
(22) 



We omitted the term (cu x , cv x ) T on the right hand sides of all systems. It is a Lie 
symmetry, corresponding to the invariance of our classification problem with respect to 
the shift of x. 



2 Discussion. 

1. Admissible transformations. Some equations in the list contain arbitrary constants 
a, P, 7, 5. Not all of them are essential. 

Let us consider, for instance, the equations (15). It is easy to see that those constants 
a and j3 which are not equal to zero can be reduced to 1 via scalings of t, x, u and v. 
In this way, (15) actually describes three different equations without parameters. These 
equations correspond toa = /3 = l, a = (3 = and a — 1, (5 — 0. 

The following parameters: a in (12) and (13), 7 in (14), both a and (3 in (17) are not 
essential in the same sense. For (14) and (16) the essential parameter is the ratio of a and 
(5. Note, that if a — (3 — then (14) coincides with the nonlinear Schrodinger equation, 
which is not a separate equation in our list. 
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For some equations from the list, there exist admissible transformations of the form 



u — > p(x, t)u + q(x, t), v — > r(x, t)v + s(x, t). (23) 

"Admissible" means that the resulting equation does not depend explicitly on x and t 
and has the same form (4). Using such admissible transformations, one can remove some 
of constants in the equations of the list. 

In particular, with the help of the transformation u — > exp(ct)w, v — > exp(— ct)v one 
can remove the terms (cm, — cv ) T in (14), (16), (17). 

Using the transformations u — > u + \t + fix, v — > v — \t — fix it is possible to remove 
(5 and 7 in (12), (13). The equation (21) can be reduced to the form 

f ut = Uxx + 6(« + v)v x - 6(« + vf + c^, , 24 , 
I v t = -v xx + 6(u + v)u x + Q(u + v) 3 + cu x 

by such a transformation and by shifts of u and v. It seems to us that the essential 
constant c was missed in the classification result of [3] . 

More general transformations are described in [4], [5] which reduce some of the equa- 
tions (8)-(22) to others in this list. For simplicity in applying our results, we will not 
rely on these non-trivial transformations, and will instead operate with the complete list 
(8)-(22). 

2. Three groups of equations. All equations of the list can be divided into three 
groups. The first group contains the so called NLS-type equations (8), (10), (13), (14), (15). 
Besides a higher symmetry (5) every such equation possesses a symmetry of the form 

{u T u xxx ~\~ ^{u, v, u X) v X) u xx , v xx ^.) (25) 
V xxx + ?p{u, V, U X) V x , U XX) V xx ). 

This is typical for equations having the Lax representations in si (2). 

The equations of the Boussinesq type form the second group (12), (20), (21), (22). They 
have no symmetries of third order. This indicates the existence of a Lax representation 
in s/(3). We have chosen the existence of the symmetry (5) as a criterion of integrability 
for (4) since the choice of the simplest ansatz (25) leads to the loss of all equations of the 
second group. 

The last group (9), (11), (16), (17), (18), (19) consists of " linearizable" equations, which 
have no higher conservation laws. Some of them seem to be new. 

Equations (16), (17) : In [18], [6] one can find a linearization procedure for (16). 
Namely, a non-local substitution 

U = -uexp (a J uvdx^j, V = vexp ^— (3 J uvdx^j (26) 

reduces it to the linear equation U t = U xx + 7?/, V t = —V xx — 7U. 

Consider equation (17). It is easy to see that it has the following symmetry 

\u T = u xxx + 3uvu xx + 6uu x v x + 3vu 2 x + 3u 2 v 2 u x , ^ 
\v T = v xxx + 3uvv xx + 6vu x v x + 3uv 2 x + 3u 2 v 2 v x , 
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for any a, (3, 7. Under a reduction u = v, (27) coincides with the well known equation (see 
[13],[11], [6]) 



u T u xxx -\- 3it *u xx + 9uu x + 3u u x , 



.2. 



which can be linearized by the substitution U = u exp J u dx. It is easy to verify that the 
following very similar substitution (cf. also with (26)) 

U = u exp ^ uv dx^j , V — v exp (^J uv dx^j (28) 

reduces (27) to U T = U xxx , V T = V xxx . The same substitution reduces (17) to a linear 
system 

U t = U xx + 2aV x + >yU, V t = -V xx - 2/3U x - 7 V. 

Generalizing the formula (28) one can find the following vector generalizations of the 
systems (17) and (27): 

' ut = u xx + 2 < u, v > u x + 2 < u, v x > u+ < u, v > 2 u+ 

2av x + /3<u,u>u + 2a<u,v>v — a < v,v > u + ju, 

v t = —v xx — 2 < u, v > v x — 2 < v , u x > v— < u, v > 2 v— 

2f3u x — a < v,v > v — 2j3 < u,v > u + j3 < u,u > v — 71;, 

Ut = u xxx + 3 < u, v > u xx + 3u < u x , v x > +3 < u, v > x u x + 3 < u, v > 2 u x , 
v T = v xxx + 3 < u, v > v xx + 3v < u x , v x > +3 < u, v > x v x + 3 < u, v > 2 v x , 

where u and v are N- dimensional vectors and < > is a scalar product. Both of them 
can be linearized just as in the scalar case: 

U = wexp J < u,v > dx, V = vexp J < u,v > dx. (29) 

In contrast with (16), (17), equations (9) and (11) are related to linear systems 

U t = U xx + Cl U x + c 2 V x + c 3 U, V t = -V xx - hV x - k 2 U x - k 3 V (30) 

via local differential substitutions. 

Equation (9): For a = [3 = 7 = the substitution is given by (3). In the general case 
such a substitution is defined by 

U x ( Cl -k 2 )U + c 2 V = V x ( Cl -c 2 )V + k 2 U 

U (U + V) 2(U + V) ' V (U + V) 2(U + V) ' 

where the constants in (30) and (9) satisfy the following conditions 

, , c 1 (c 1 -c 2 -k 2 ) 
ki = ci, fc 3 = c 3 = , 
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2a = —k 2 (ci + c 2 ), 2f3 = — c 2 (ci + fc 2 ), 27 = ci/c 2 c 2 . 
Equation (11): The substitution is of the form 

c 2 U x c\ a 

u = — +r v = v x /v + -. 

The relations between constants are the following 

h = ci, k 2 = c 2 , k 3 = 0, 

2a = -c\ - 2c 2 + 2c 3 , /3 = -cic 2 , 47 = c 2 (2ci + c 2 - 2c 3 ). 

Equation (18), (19): These equations are related to triangular systems. Equations 
(19) have been obtained in [17]. The substitution 

U x V 

U= w^ V= VU ] 

found first by Marihin [15], links equations (19) with 

U t = U xx + 2V 2 , V t = -V xx + aV. 

The last system is linear in the following sense. To find V we need to solve a linear 
equation. For a given function V, the function U satisfies a linear equation with variable 
coefficients. 

Similarly the following substitution 

u = -U- 2/ *V- l U X) v = U- 1/3 V, 
3 

reduces (18) to 

U t = U xx - 2V- 1 V X U X + 3aU + 3V 3 , V t = -V xx . 

3. Master-symmetries. For all equations (8)-(22) we have found all symmetries of 
order less or equal than four using the computer program LiePde [20]. It turns out that 
many of the equations have symmetries depending on t and x explicitly. For example, 
computing for equation (16) the general symmetry of n order of the form 

u T = P{t)u nx + . . . , v T = -P{t)v nx + ... 

the polynomial P(t) is an arbitrary polynomial of degree n. 

It is well known that symmetries which are linearly depending on t and x are closely 
related to local master-symmetries [14]. The equations (8)-(ll), (15)-(17) have such 
symmetries. To obtain the master-symmetries one has to simply put t equal to zero in 
these time-dependent symmetries. 
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The resulting master-symmetry is of the form 

f u T = 2x (u xx + F(u, v, u x , v x )) + f(u, v, u x , v x ), , . 

\v T = 2x (-v xx +G(u,v,u x ,v x )) +g(u,v, u x , v x ), 

where F and G are the right hand sides of the corresponding equation (1) and / and g 
are given by the following list 



(8) 


/ = 


u 2 + 3uv + Au x , 


9 = 


v 2 + 3uv — Av x , 


(9) 


/ = 


2u 2 + Auv + (3 + 3u x , 


9 = 


—2v 2 — Auv — a — 3v x , 


(10) 


/ = 


iuv + 5u x , 


9 = 


v 2 + 4u — 5v x , 


(11) 


/ = 


Aaiv + P + 3u x , 


9 = 


—2v 2 — 2u — a — 3v x , 


(15) 


/ = 


2u 2 v + 2av + 3u x , 


9 = 


2uv 2 + 2(3u - 3v x , 


(16) 


/ = 


2au 2 v — 2jxu + 2u x , 


9 = 


2(3v 2 u + 2>yxv-2v x , 


(17) 


/ = 


2u 2 v + 2av + 2u x , 


9 = 


-2uv 2 - 2f3u - 2v x . 
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